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Abstract 

This paper describes an efficient and unique method for computing the shortwave 
direct radiative effect (DRE) of aerosol residing above low-level liquid-phase clouds 
using CALIOP and MODIS data. It addresses the overlap of aerosol and cloud rigorously 
by utilizing the joint histogram of cloud optical depth and cloud top pressure while also 
accounting for subgrid-scale variations of aerosols. The method is computationally 
efficient because of its use of grid-level cloud and aerosol statistics, instead of pixel-level 
products, and a pre-computed look-up table based on radiative transfer calculations. We 
verify that for smoke over the southeast Atlantic Ocean the method yields a seasonal 
mean instantaneous (approximately 1:30PM local time) shortwave DRE of above cloud 
aerosol (ACA) that generally agrees with more rigorous pixel-level computation within 
4%. We also estimate the impact of potential CALIOP aerosol optical depth (AOD) 
retrieval bias of ACA on DRE. We find that the regional and seasonal mean 
instantaneous DRE of ACA over southeast Atlantic Ocean would increase, from the 
original value of 6.4 W m" 2 based on operational CALIOP AOD to 9.6 W m 2 if CALIOP 
AOD retrieval are biased low by a factor of 1.5 (Meyer et ah, 2013) and further to 30.9 W 
m 2 if CALIOP AOD retrieval are biased low by a factor of 5 as suggested in (Jethva et 
ah, 2014). In contrast, the instantaneous ACA radiative forcing efficiency (RFE) remains 
relatively invariant in all cases at about 53 W nf 2 AOD' 1 , suggesting a near linear relation 
between the instantaneous RFE and AOD. We also compute the annual mean 
instantaneous shortwave DRE of light-absorbing aerosols (i.e., smoke and polluted dust) 
over global oceans based on 4 years of CALIOP and MODIS data. We find that the 
variability of the annual mean shortwave DRE of above-cloud light-absorbing aerosol is 
mainly driven by the optical depth of the underlying clouds. While we demonstrate our 
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1 method using CALIOP and MODIS data, it can also be extended to other satellite data 

2 sets, as well as climate model outputs. 
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1. Introduction 


The shortwave direct radiative effect (DRE) of aerosols at the top of the atmosphere 
(TOA) is strongly dependent on the reflectance of the underlying surface. Over dark 
surfaces (e.g. ocean, vegetated land), the scattering effect of aerosols is generally 
dominant, leading to negative DRE (i.e., cooling) at TOA (Yu et al., 2006). In contrast, 
when light-absorbing aerosols occur above clouds or other bright surfaces (such as snow, 
ice, desert), aerosol absorption is significantly amplified by cloud or surface reflection, 
offsetting or even exceeding the scattering effect of the aerosol, leading to a less negative 
or even positive (i.e., wanning) TOA DRE (Abel et al., 2005; Keil and Haywood, 2003; 
Twomey, 1977). Therefore, in order to understand the full complexity of aerosol radiative 
effects on climate, it is important to quantify the DRE under both clear-sky and cloudy- 
sky conditions. Although the DRE of aerosols in clear-sky regions has been extensively 
studied and is relatively well constrained based on advanced satellite remote sensing 
measurements acquired in the last decade (Yu et al., 2006), the cloudy-sky DRE is 
generally assumed to be negligible or simulated by models (Schulz et al., 2006). 

Currently model simulations shows a large inter-model spread in cloudy-sky DRE 
(Schulz et al., 2006), which results from inter-model differences in both aerosol and 
cloud properties (Schulz et al., 2006; Stier et al., 2013). Therefore, there is a clear need 
for an observational constraint on the DRE of above-cloud aerosol (ACA). 

Recent advances in satellite remote sensing techniques have provided an 
unprecedented opportunity for studying the DRE of ACA. In particular, the availability of 
measurements from the space-borne Cloud-Aerosol Lidar with Orthogonal Polarization 
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(CALIOP) sensor onboard NASA’s Cloud-Aerosol Lidar and Infrared Pathfinder 
Satellite Observations (CALIPSO) satellite has provided a revolutionary global view of 
the vertical distribution of aerosols and clouds (e.g., Winker et ah, 2013). Using CALIOP 
aerosol and cloud layer products, Devasthale and Thomas (2011) found frequent 
occurrences of aerosols residing above low-level clouds in several regions of the globe. 

In particular, they found a high frequency of smoke occurrence over low clouds in the 
southeast Atlantic, western coasts of South America (e.g., Columbia, Ecuador, and Peru) 
and southern Asia. These authors also found a high frequency of natural and polluted dust 
aerosols overlapping low clouds off the western coasts of Saharan Africa in boreal 
summer and over boundary layer clouds in the eastern coast of China in boreal spring 
(see Fig. 3 of Devasthale and Thomas, 2011). 

CALIOP measurements of ACA properties, in combination with satellite cloud 
products from, for example, the Moderate Resolution Imaging Spectroradiometer 
(MODIS), have been used in several recent studies to derive the DRE of ACA with 
radiative transfer simulations (e.g., Chand et ah, 2009; Costantino and Breon, 2013b; 
Meyer et ah, 2013; Oikawa et ah, 2013). (Chand et al., 2009) used CALIOP above-cloud 
AOD retrievals (Chand et al., 2008) and Terra-MODIS cloud products, both aggregated 
to 5° gridded monthly means, to calculate the radiative effects of smoke transported 
above the low-level stratocumulus deck in the southeastern Atlantic. A major point made 
in this study was that the all-sky DRE of elevated light-absorbing aerosols, such as 
transported smoke, is strongly modulated by the underlying cloud properties. However, 
the spatial-temporal aggregation of both CALIOP and MODIS data to coarse gridded 
monthly means obscures the potential influence of cloud and aerosol variability on the 
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DRE. In particular, using grid box mean cloud optical depth for DRE calculation might 
lead to biases in DRE due to the plane-parallel albedo bias (Oreopoulos et ah, 2007). 
Moreover, the MODIS level-3 aggregation algorithm samples all liquid water clouds, 
regardless of possible retrieval contamination by ACA. As a result, the total population of 
liquid water clouds in the MODIS level-3 products (daily or monthly) may be 
significantly different from that of below-aerosol-only cloud population. Therefore, using 
level-3 MODIS products without distinguishing below-aerosol-only from total cloud 
population can potentially lead to significant errors. The problem could be further 
complicated by biases in MODIS cloud retrievals associated with the presence of 
overlying light-absorbing aerosols. When a cloud-pixel is contaminated by overlying 
light-absorbing aerosols the MODIS cloud optical depth (COD) retrieval is generally 
biased low (e.g., Coddington et ah, 2010; Haywood et ah, 2004; Jethva et ah, 2013; 
Wilcox, 2010), an effect not considered in most previous studies (e.g., Chand et ah, 2009; 
Costantino and Breon, 2013b; Oikawa et ah, 2013). Most recently however, (Meyer et ah, 
2013) collocated CALIOP above-cloud AOD and Aqua-MODIS cloud properties at the 
pixel level, and the DRE was then computed at these individual collocated pixels. They 
found that correcting the MODIS COD bias due to ACA contamination can lead to a 
more positive ACA DRE. Such rigorous collocation has obvious advantages as it takes 
into account the sub-grid variability of clouds and aerosols, but is on the other hand 
computationally expensive since it requires large amounts of pixel-level data that make 
global scale and multiyear studies challenging. 

The objective of this paper is to describe a novel method for computing the DRE of 
ACA. This method attempts to balance the need for computational efficiency with the 
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need for rigorous treatment of aerosol-cloud overlap and small-scale variability of aerosol 
and clouds. Our method has several unique features: 1) it takes sub-grid scale cloud and 
aerosol variation into account in DRE computations; 2) it treats the overlap of aerosol and 
cloud rigorously by utilizing the joint histogram of COD and cloud top pressure (CTP) in 
the MODIS level-3 product; 3) it is computationally efficient because of the use of a pre- 
computed look-up table of ACA DRE. 

In the following sections, we briefly introduce the CALIOP and MODIS data used 
(Section 2), describe the key assumptions and features of the novel method (Section 3), 
validate it through comparison with pixel-level computations as in (Meyer et ah, 2013) 
(Section 4), and conclude with a summary and discussion (Section 5). 

2. Satellite Data 

In (Meyer et al., 2013), the MODIS level-2 cloud product is collocated with CALIOP 
level-2 aerosol product for every pixel along the CALIOP track and the computation of 
instantaneous DRE is performed pixel-by-pixel. Then, the pixel-level DRE results are 
aggregated on a latitude-longitude grid for climatological study. If only the grid-level 
DRE is of interest, the pixel-by-pixel computation of DRE may not efficient because of 
redundant computations. For example, if two pixels with the same above-cloud AOD and 
below-cloud COD occur within the same grid-box, they evidently have the same ACA 
DRE, but the radiative transfer computation is nevertheless performed twice in the pixel- 
by-pixel method. As shown in Section 3, a more efficient way is to compute the DRE 
statistically using the probability density function (PDF) of above-cloud AOD and below- 
cloud COD. In this study, we use the CALIOP level-2 aerosol and cloud layer product 
(V3.01) to derive the statistics of ACA properties and the MODIS level-3 daily cloud 
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product for cloud property statistics. It important to note that our method is not limited to 
CALIOP and MODIS products, but also applicable to other satellite data sets, such as 
above-cloud aerosol retrievals from POLDER (POLarization and Directionality of 
the Earth's Reflectances) (Waquet et al, 2009) and OMI (Ozone Monitoring 
Instrument)(Torres et al., 2012), and cloud retrievals from ISCCP (International Satellite 
Cloud Climatology Project) (Rossow and Schiffer, 1999) and SEVIRI (Spinning 
Enhanced Visible and Infrared Imager)(Schulz et al., 2009), as well as outputs from 
general circulation models and chemical transport models. 

2.1. CALIOP level-2 aerosol and cloud layer products 

Since its launch in 2006, the space-borne lidar CALIOP has continuously acquired, 
with near global (albeit instantaneously sparse) coverage, attenuated backscatter 
measurements at 532 nm and 1064 nm, including linear depolarization information at 
532nm (Winker et al., 2009). The CALIOP level-2 retrieval algorithm consists of several 
steps. First, a “feature finder” algorithm and cloud-aerosol discrimination (CAD) 
algorithm are used to detect aerosol and cloud layers, and record their top and bottom 
heights and layer integrated properties (Vaughan et al., 2009). Second, the detected 
aerosol layers are further classified into six sub-types (i.e., polluted continental, biomass 
burning, desert dust, polluted dust, clean continental and marine) (Omar et al., 2009) and 
the detected cloud layers are assigned different thennodynamic phases (Hu et al., 2007a) 
based on the observed backscatter, color ratio and depolarization ratio. Third, a priori 
lidar ratios, pre-selected based on aerosol sub-type and cloud phase, are used to derive the 
extinction of an aerosol or cloud layer from the attenuated backscatter profile (Young and 
Vaughan, 2008). 
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In this study, we use CALIOP level-2 version 3.01 aerosol and cloud layer products at 
a nominal 5 km horizontal resolution (i.e., CAL_LID_L2_05kmALay and 
CAL_LID_L2_05kmCLay) for aerosol-cloud overlap detection, and for information on 
aerosol layer properties, including type, aerosol optical depth (AOD), and layer top and 
bottom height. In addition to physical properties, the CALIOP layer products also provide 
various metrics and flags on data quality assurance. These include CAD score (Liu et ah, 
2009), horizontal averaging scale, extinction quality control (QC) flag, and estimated 
uncertainty of layer AOD. In this study, we apply these metrics following best practices 
provided by the CALIPSO science team to screen for reliable retrievals (e.g., Winker et 
al., 2013) (see Table 1). 

It should be noted here that the current version of CALIOP operational aerosol 
retrieval algorithm (V3.01) appears to significantly underestimate the AOD of above- 
cloud aerosol layer according to recent studies (Jethva et al., 2014; Liu et al., 2013; 
Waquet et al., 2013b). The main reason is that after strong attenuation by the upper part 
of an aerosol layer, the 532 nm attenuated backscatter of the lower part of aerosol layer is 
often too small. As a result, it is not included in the current CALIOP retrieval, leading to 
AOD that is biased low (Jethva et al., 2014; Liu et al., 2013). At the moment, the 
CALIOP operational team is investigating the possibility of using the algorithm described 
in (Chand et al., 2008; Hu et al., 2007b) for ACA retrievals (Liu et al., 2013). This 
alternate method utilizes the reflected lidar signal from the cloud layer underneath to 
derive the two-way transmittance and thereby the AOD of the ACA layer. Because the 
backscatter of a cloud layer is usually very strong, the two-way transmittance method is 
less affected by the strong attenuation of the ACA layer and is therefore expected to 
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alleviate the aforementioned problem. Lidar based AOD retrievals are also known to 
suffer from other issues, such as the background solar noise during daytime. These issues 
are beyond the scope of this study, but are neverthelss discussed in the uncertainty 
analysis of Section 3.4. 

In addition to retrieval errors and uncertainties, another limitation of CALIOP data is 
the small sampling rate (i.e., only along track). In order to compute the DRE of ACA 
over a given latitude-longitude grid box, we assume that the aerosol property statistics 
retrieved by CALIOP along its narrow track represent the statistics over the whole grid 
box, i.e., that AOD PDFs are identical. This assumption constitutes an uncertainty in our 
DRE computation. With no other data available to compare, it is difficult to determine the 
size of this uncertainty. Recently, several novel methods have been developed to retrieve 
ACA properties from passive sensor observations (Jethva et ah, 2013; Torres et ah, 2012; 
Waquet et ah, 2009; 2013a), which will help improve our understanding of the sub-grid 
ACA variability when they become available to public. 

Finally, we emphasize two more points. First, none of the aforementioned problems 
with CALIOP data, e.g., smoke AOD bias, retrieval uncertainties, and small sampling 
rate, are unique to our method. Any method that uses CALIOP data faces the same 
challenges. Second, our method is not limited only to CALIOP data. We choose to use 
CALIOP product in this study solely because it is the only publically available ACA 
product at the moment. Our method can also be applied to other ACA retrieval products, 
based on for example, POLDER (Waquet et al., 2009), MODIS (Jethva et al., 2013), 

OMI (Torres et al., 2012) observations when they become available to public. In fact, as 
discussed later, the advantage of our method in terms of computational efficiency is even 
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greater when applied on retrievals from passive sensors. 


2.2. MODIS daily level-3 cloud property product 

This study computes the grid-level ACA DRE using the statistics of aerosol and cloud 
properties, instead of pixel-by-pixel computation as in (Meyer et ah, 2013). We use the 
Collection 5 (C5) Aqua MODIS level-3 Daily gridded Atmosphere product MYD08 D3 
for the statistics of cloud properties and other parameters, such as solar zenith angle, 
needed for ACA DRE computations. 

The MODIS level-3 (i.e., grid-level) product contains statistics computed from a set 
of level-2 (i.e., pixel-level) MODIS granules. As summarized in (Platnick et ah, 2003), 
the operational level-2 MODIS cloud product provides cloud masking (Ackennan et al., 
1998), cloud top height retrieval based on CO 2 slicing or the infrared window method 
(Menzel et al., 1983), cloud top thennodynamic phase determination (Menzel et al., 

2006), and cloud optical and microphysical property retrieval based on the bi-spectral 
solar reflectance method (Nakajima and King, 1990). In addition to these cloud 
parameters, the level-2 products also provide pixel-level runtime Quality Assessment 
(QA) information, which includes product quality as well as processing path information. 
All MODIS level-2 atmosphere products, including the cloud, aerosol and water vapor 
products, are aggregated to 1° spatial resolution on a daily (product name MYD08 D3 
for Aqua MODIS), eight-day (MYD08 E3), and monthly (MYD08 E3) basis. 
Aggregations include a variety of scalar statistical information (mean, standard deviation, 
max/min occurrences) and histograms (marginal and joint). A particularly useful level-3 
cloud product for this study is the daily joint histogram of COD vs. CTP, derived using 
daily counts of successful daytime level-2 pixel retrievals that fall into each joint COD- 


10 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 


CTP bin. Eleven COD bins, ranging from 0 to 100, and 13 CTP bins, ranging from 200 to 
1000 mb, comprise the histogram. As discussed below, the COD-CTP joint histogram 
allows for identification of the portion of the cloud population that lies beneath the 
aerosol layer found by CALIOP, as well as the corresponding COD probability 
distribution needed for DRE estimation. In addition to the COD-CTP joint histogram, we 
also use the gridded mean solar and sensor zenith angles for calculating DRE and 
correcting the COD bias due to the presence of ACA. 

It should be noted that the level-3 daily product MYD08 D3 contains statistics 
computed from a set of level-2 MODIS granules that theoretically span a 24-hour interval 
(Hubanks et ah, 2008). However, for cloud parameters retrieved only during daytime, 
such as COD and cloud droplet effective radius (CER), only daytime level-2 files are 
used to compute the level-3 daily statistics. These are called daytime only SDSs 
(Scientific Data Sets) in level-3 products. Strictly speaking, the daytime only SDSs of 
only those 1° gridcells between approximately 23° N and 23° S come from a single 
MODIS overpass. The tropical South East Atlantic region, where transported smoke 
aerosols are often observed above low-level stratocumulus clouds, is within this range 
(about 10° N~30° S see Figure 3). The COD statistics in MYD08 D3 product for this 
region are therefore derived from a single Aqua-MODIS overpass that can be collocated 
with CALIOP observations (see Section 3.1 for details on collocation). The DRE 
computed based on the collocated dataset is therefore instantaneous DRE at Aqua 
crossing time (1:30PM) that are directly comparable to the pixel-by-pixel results 
in(Meyer et ah, 2013). Poleward of 23°, MYD08_D3 statistics are derived from 
averaging several overlapping orbits approximately 100 minutes apart(Hubanks et ah, 
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2008). As a result, strictly speaking for mid and high latitude regions the DRE computed 
based on MYD08 D3 data is not instantaneous DRE. We emphasize that this is not a 
limitation of our method, but an inherent characteristic of the MODIS level-3 product. 

3. Methodology 

3.1. Theoretical basis 

As in previous investigations (e.g., Chand et ah, 2008; 2009; Costantino and Breon, 
2013b; Meyer et ah, 2013), we focus on the simplest case of overlapping aerosol and 
cloud, i.e., a single layer of aerosol overlying a single layer of low-level liquid-phase 
clouds, which is commonly observed in many regions of the globe (Devasthale and 
Thomas, 2011). More complex situations cerianly exist, such as an aerosol layer located 
in between high and low cloud, or an aerosol layer overlying multiple layers of clouds. 
However, identification of such situations are either beyond the detection capabilities of 
CALIOP or relatively rare (Devasthale and Thomas, 201 1). As such, they are not 
considered here and left for future research. 

To illustrate the theoretical foundation of the method, consider the schematic example 
in Figure 1. For a given grid box (e.g., l°xl° in case of MODIS level-3 data), the gridded 
mean instantaneous broadband shortwave DRE ( ^ DRE ^ ACA ) averaged over all ACA 
pixels within the grid box is given by: 

(DRE) Aa = ff DRE{1„1 Mi (1) 

where p( T c ,T a ) is the joint probability density function (PDF) of the above-cloud AOD 
at 532 nm ( T ) and below-aerosol COD ( T c ) of ACA pixels. We note that, in addition to 
T a , DRE also depends on the the spectral variation of aerosol and cloud optical depth, 
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spectral single scattering albedo and asymmetry factor, wavelength dependencies not 
explicitly shown in this equation. These properties are computed using a Mie scattering 
code (Wiscombe, 1980) based on the aerosol model described in (Meyer et ah, 2013). 

The dependencies on solar zenith angle, surface reflectance, cloud particle effective 
radius, and atmospheric profile are also omitted from the equation; solar zenith angle and 
surface reflectance are expected to have only minor variation within the grid box, while 
the impact of cloud particle effective radius and atmospheric profile on shortwave DRE is 
relatively small. Since p ( r ,T ) describes the covariation of aerosols and clouds for the 

ACA pixels, it should ideally be derived from collocated CALIOP aerosol and MODIS 
cloud retrievals at pixel level as in (Meyer et ah, 2013). This requires large amounts of 
pixel-level data, however, as one month of global daytime C5 MODIS level-2 cloud 
products in HDF format are roughly 150 Gigabytes. Therefore, pixel-level collocation 
and radiative transfer simulation are too computationally expensive and cumbersome for 
multiyear global studies. 

A key assumption in our method, which allows us to avoid tedious pixel-level 
collocation, is that the sub-grid level instantaneous spatial distribution of above-cloud 
AOD is statistically independent from the sub-grid level instantaneous spatial distribution 
of below-aerosol COD. Under this assumption, p( T c ,T a ) = p(l c )- p(T a ) and Eq. (1) 
reduces to: 


DRE )aCA = \l\^ DRE ^c^a)P^c) dT c 


( 2 ) 


where p{ t c ) and p{ r a ) are the PDF of instantaneous below-aerosol COD T . and above- 


cloud AOD T a , respectively, of ACA pixels. The advantage of Eq. (2) is that it allows 
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p{ r ,) and p{ T u ) to be derived separately and independently. This assumption is 

reasonable considering that transported ACAs and low-level boundary layer clouds are 
usually well separated vertically (Devasthale and Thomas, 2011) and controlled by 
different meteorological conditions. The potential coupling between the two is that 
overlying absorbing aerosols could influence the evolution of clouds through changing 
atmospheric stratification (Wilcox, 2010). However, a recent observational study 
(Costantino and Breon, 2013a) found no correlation between above-cloud AOD and 
below-aerosol COD, although correlations are found between AOD and cloud droplet 
effect radius, as well as liquid water path. Moreover, it is important to stress that our 
assumption is that the instantaneous above-cloud AOD and below-aerosol COD are 
independent at sub-grid scale. This assumption does not rule out the possibility that AOD 
and COD could be correlated at longer temporal (e.g., seasonal) and/or larger spatial (e.g., 
regional) scale through the thennodynamic and radiative coupling (Wilcox, 2010; 2012). 
Finally, as shown in section 4, when we compare the DRE derived from pixel-level 
collocation (i.e., based on Eq. (1)) with that from independent sampling of p(x c ) and p(x a ) 
(i.e., based on Eq. (2)) the agreement is very good. 

In our method, the PDF of above-cloud AOD p ( r ) is derived from the CALIOP 

5km aerosol and cloud layer products through the following steps: 1) for each 5km 
CALIOP profile that falls within a given latitude-longitude grid box, we first search for 
an aerosol layer; 2) if an aerosol layer is detected and the quality metrics pass the quality 
assurance criteria summarized in Table 1, we then proceed to check for the presence of an 
underlying liquid-phase cloud layer within the profile using the CALIOP cloud layer 
product; 3) if a cloud layer is present, the AOD of the aerosol layer is recorded for the 
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derivation of the p( r ) of the grid box. The bottom height of the aerosol layer is also 

recorded to derive the grid mean aerosol layer bottom height. Once all of the CALIOP 
profiles within the grid box are processed, we obtain the PDF of the above-cloud 
AOD p(T a ) and the mean aerosol layer bottom pressure ( P ho!lom ) ■ 

As schematically illustrated in Figure 1, the PDF of below-aerosol COD p { r ) is 

derived from the joint histogram of cloud optical depth and cloud top pressure (COD- 
CTP joint histogram) in the MODIS daily level-3 product, using the grid mean aerosol 

layer bottom pressure ( P holtom ) derived above. For a given grid box, we first identify the 
population of liquid-phase clouds below the pressure level ( P hotlom ) ■ This subset, together 
with the AOD PDF p ( t ) , is then used to calculate DRE according to Eq. (2). 

In this study, we focus on the computation of instantaneous DRE. To obtain diurnally 
averaged DRE, technically speaking one would simply need to integrate over time the 
instantaneous DRE. However, it is important to note that, in addition to diurnal variation 
of solar zenith angle, aerosol and cloud properties may also have significant diurnal 
cycles. In fact, it is known that the low cloud fraction over stratocumulus regimes, such 
as the South East Atlantic region, have a strong diurnal cycle (15-35% of diurnal mean 
value) driven by cloud solar absorption(Wood et ah, 2002). A recent study by (Min and 
Zhang, n.d.) indicates that using a constant cloud fraction based on Aqua-MODIS 
observations tends to result in significantly underestimated diurnal mean DRE even if the 
diurnal variation of solar zenith angle is considered in the computation. Therefore, the 
challenge for diurnal DRE computation is to obtain high-temporal-resolution cloud and 
aerosol observations. In this regard, the SEVIRI on the geostationary satellite MSG 
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(Meteosat Second Generation) may be ideal for a diurnal cycle study since it provides 
aerosol and cloud observations over the South East Atlantic region at 15 minute intervals. 
However, a pixel-by-pixel computational approach would be too expensive to be used for 
the diurnal DRE computation based on high-temporal (15 minutes) and high-spatial (4 
km at nadir) resolution of SEVERI data. Our method, on the other hand, meets well the 
demand for computational efficiency. 

3.2. DRE Look-up Tables 

To speed up calculations, we use pre-computed aerosol-type specific look-up-tables 
(LUTs), instead of online radiative transfer computation, when deriving the DRE of AC A. 
The concept of our LUTs is somewhat similar to the “radiative kernels” described in 
(Hartmann et ah, 2001) and (Zelinka et ah, 2012) for computing cloud radiative 
feedbacks. The LUT for each aerosol type consists of DREs at both TOA and surface (not 
used in this study) for various combinations of AOD, COD, CTP and solar zenith 
conditions. As such, once the aerosol type and AOD are known from CALIOP and COD, 
CTP and solar zenith angle are known from MODIS, the corresponding DRE can be 
obtained through LUT interpolation. Note that the CALIOP only provides AOD at lidar 
wavelengths (e.g., 532 nm and 1064 nm) for each aerosol type. Therefore, radiative 
transfer model-appropriate narrowband aerosol scattering properties, namely AOD, 
single-scattering albedo and asymmetry factor, are needed for the development of the 
LUT. The current version of LUT focuses on light-absorbing aerosols (e.g., smoke and 
polluted dust). In order to validate our method with more rigorous pixel-level 
computations, we adopt the narrowband aerosol optical properties of (Meyer et ah, 2013), 
who used the same radiative transfer code, in the computation of the current LUT. The 
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aerosol model in (Meyer et al., 2013) is extended from an absorbing aerosol model 
developed for the MODIS Collection 5 Aerosol Product (MOD04) (see Table 4 of Levy 
et al., 2009). The MOD04 aerosol models define aerosol size distributions and refractive 


indices based solely on prescribed AOD at 550 nm (MODIS band 4; note that the 
absorbing aerosol model used here assumes a constant index of refraction, 1 .5 1— 0.02i, at 
all wavelengths). At AOD=0.5 (550 nm), the single-scattering albedo of this model is 
about 0.9 over the visible spectral region (see Figure 7 of Meyer et al., 2013), which is in 
the range of previously reported values (e.g., Keil and Haywood, 2003; Myhre et al., 
2003). The current AOD bins (at 550 nm) in the LUT range from 0.05 to 1.5, which 
covers most of the above-cloud AOD observed by CALIOP. The current COD bins, 
logarithmically spaced, range from 0.1 to 300. Following the MODIS level-3 data, the 
thirteen CTP bins range from lOOOmb to 200mb. The solar zenith angle bins range from 0 
to 80 degree. Radiative transfer computations are carried out using the RRTM-SW model 
(Clough et al., 2005; Iacono et al., 2008). Lambertian ocean surface reflectance is set to 
5%. Cloud droplet effective radius is fixed at 15 pm, which is close to the global mean 
value over oceans observed by MODIS (King et al., 2013). This value of effective radius 
is also used to convert the MODIS visible COD to liquid water path used as input to 
RRTM-SW. Liquid cloud optical properties are calculated internally by RRTM. For 
atmospheric profiles of water vapor and temperature, we use NCEP R1 reanalysis data 
(Kistler et al., 2001) averaged both zonally and annually. Our sensitivity tests indicate 


that the shortwave DRE of ACA is largely insensitive to cloud effective radiu; 


or 


atmospheric profiles. 
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3.3. Cloud Optical Depth Correction 

As noted in previous studies (Coddington et al., 2010; e.g., Haywood et al., 2004), 
when a cloudy MODIS pixel is contaminated by overlying light-absorbing aerosols the 
COD retrieval is generally biased low. We have developed a fast COD correction scheme 
to account for the COD retrieval bias due to ACA in our DRE computation, which is 
illustrated in Figure 2. This scheme requires both the cloud reflectance LUT for clouds 
without ACA, for which we use the MODIS operational LUT, and clouds with ACA, for 
which we use the one developed by Meyer et al. (2013). In the operational MODIS 
retrieval, the reflectance LUT of cloud without ACA is used to interpret the reflectance of 
all clouds, including those affected by ACA. Based on this fact, we first infer the 
“observed” cloud reflectance (after atmospheric correction) by interpolating the 
reflectance LUT of cloud without ACA corresponding to the biased COD. Then, we use 
the “observed” cloud reflectance and ACA-affected LUT (derived based on CALIOP 
AOD) to detennine the corrected COD. This COD correction process is performed for 
every combination of COD bin in p(x c ) and AOD bin in p(x a ). In the final step we 
resample the corrected CODs to obtain the corrected p(x c ). 

It should be noted that because different aerosol type may have different impact 
on MODIS COD retrievals, the above COD correction process is aerosol-type dependent. 
In this study, we use light-absorbing aerosols as example to illustrate our method and for 
validation purposes we use the aerosol model developed by (Meyer et al., 2013) for the 
development of LUTs for DRE computation and COD correction. However, the LUTs 
can be easily extended to other aerosol models. In fact, as part of ongoing research, we 
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are extending our LUTs to include all six operational CALIOP aerosol models as 
described in (Omar et al, 2009). 

3.4. Uncertainty Analysis 

Several recent studies suggest that the current operational CALIOP product tends to 
underestimate the above-cloud AOD. Meyer et al. (2013) found that the daytime 
CALIOP AOD retrievals are systematically smaller than the nighttime retrievals, 
probably due to the daytime solar background issue. In the light of this finding, Meyer et 
al. (2013) increased the CALIOP AOD retrievals by a factor of 1.5 to account for the 
impact of potential AOD bias on DRE of ACA. A more recent case study by (Jethva et 
al., 2014) suggests that CALIOP ACA AOD retrievals are biased low by a factor of 5 or 
even more compared with other retrievals, although the generality of this finding needs to 
be further tested with larger samples. While a rigorous analysis uncertainty analysis of 
CALIOP AOD product is beyond the scope of this study, it is nevertheless reasonable to 
assume that the current CALIOP retrievals provide a lower limit to the ACA AOD. In the 
uncertainty analysis presented in the next section, we carry out two sensitivity tests to 
estimate the potential impacts of CALIOP AOD bias on DRE computation. We multiply 
CALIOP AOD values by a factor of 1.5 in the first test following (Meyer et al., 2013) and 
by a factor of 5 in the second as suggested in (Jethva et al., 2014). 

Once the magnitude of the uncertainties in the input data is prescribed, the 
consequential impact on DRE can be easily estimated in our method as follows. First, in 
addition to the p( T. ; ) based on the original CALIOP data, we also derive the perturbed 

PDF p( t a ) by perturbing the original data according to pre-defined uncertainties (i.e., by 
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increasing the original values by a factor of 1.5 or 5). Then, the impact of input 
uncertainty on ACA DRE can be estimated by comparing the DREs computed with the 

original vs. perturbed PDF (i.e., ^ DRE 'j ACA vs. l^DRE^j )C( ). Note that l^DRE ^ ACA and 

DRE ^ ACA can be obtained in a single computation because they both represent integrals 

over DRE(t c ,t ) , only with different weights. In this regard, our method is much more 

efficient than the pixel-by-pixel method, in which uncertainty must be estimated by 
perturbing individual pixels. 

4. Implementation and validation of new DRE estimation scheme 

Each year during austral winter, dry season biomass burning activities throughout 
southern Africa inject large amounts of smoke into the troposphere (Eck et ah, 2003; 
Ichoku et ah, 2003; Myhre et ah, 2003). Prevailing easterly winds during this season 
often transport the smoke westward off the continent, over the ocean, where extensive 
marine boundary layer clouds persist for most of the year. Under the descending branch 
of the Hadley cell, the air mass above the boundary layer is quite dry. Due to the lack of 
efficient wet scavenging, the transported aerosol layers can remain suspended in the 
atmosphere for days, creating a near-persistent smoke layer above the stratocumulus deck 
over the southeastern Atlantic Ocean (Chand et ah, 2009; Devasthale and Thomas, 2011; 
Keil and Haywood, 2003; Wilcox, 2010). 

To validate our method, we have compared the DRE of above-cloud light-absorbing 

aerosols in this region with pixel-level computations from (Meyer et ah, 2013). Figure 3a 

shows the seasonal mean (August/September 2007-2011) instantaneous TOA DRE of 

above-cloud smoke and polluted dust based on the pixel-level computations from (Meyer 
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et al., 2013). Figure 3b shows the corresponding instantaneous TOA aerosol radiative 
forcing efficiency (RFE) defined as the DRE per unit AOD. The DRE and RFE results 
computed using our method described in the previous section are shown in Figure 3 c and 
Figure 3d, respectively. Evidently, both DRE and RFE computed using our new method 
agree closely with the pixel-level computations. Figure 4 shows the meridional mean 
DRE and RFE for the region using the results in Figure 3. Not surprisingly, the outcomes 
of the two methods are almost identical. Note that the CODs used in the computations for 
Figure 3 and Figure 4 are directly from the MODIS products without COD correction. 

We have also compared the DRE and RFE from the two methods using the corrected 
COD and achieved again very good agreement (not shown because of close resemblance 
to Figure 3 and Figure 4). The seasonal and regional mean DRE and RFE, based on the 
corrected COD, from the pixel-level computation method in(Meyer et al., 2013) are 6.63 
W m' 2 and 55.97 W m' 2 AOD' 1 , respectively (see Table 2). The corresponding values 
from our new method are 6.39 W m' 2 and 53.77 W m' 2 AOD' 1 , respectively. 

As previously mentioned, to estimate potential bias in CALIOP ACA AOD retrieval 
on our DRE computation, we carried out two sensitivity tests. We increased CALIOP 
AOD values by a factor of 1.5 in one test following (Meyer et al., 2013) (hereafter 
refered to as “xl.5 test”) and by a factor of 5 in another as suggested in (Jethva et al., 
2014) (“x5.0 test” hereafter). In both cases, we corrected the MODIS COD retrievals 
based on the scaled AOD. The regional and seasonal mean DRE of ACA increases, from 
the original value of 6.39 W m" 2 to 9.57 W m' 2 in the xl.5 test and to 30.87 W m' 2 in the 
x5.0 test. We have to note that this is a very rough estimate. Nevertheless, the DRE based 
on the xl.5 scaling of CALIOP AOD seems to agree reasonably with the value, 9.2±6.6 
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W m" 2 , reported in an independent study by Wilcox (Wilcox, 2012). Interestingly, the 
scaling of AOD has little impact on RFE in both cases (53.09 W m" 2 AOD' 1 in the xl.5 
case and 51.24 W irf 2 AOD' 1 in the x5.0 case), which apparently suggests a near linear 
relationship between DRE and AOD as also noted in (Meyer et ah, 2013) and (Wilcox, 
2012) (see his Figure 5). 


In summary, as shown clearly in Figure 3, Figure 4 and Table 2, the DRE inferred 
from our new method agrees very well with the pixel-level computations. Furthermore, 
the difference between the two methods is much smaller than, for example, the 
uncertainty associated with CALIOP retrieval biases. 

It is worthwhile to clarify again that the results shown in Figure 3 are seasonal mean 
instantaneous DRE at A-Train crossing time (1:30PM local time) based on CALIOP 
above-cloud AOD and corrected Aqua MODIS below-aerosol COD retrievals. Moreover, 
the aerosol model described in (Meyer et ah, 2013) is used in this study. All these factors 
should be considered when comparing the results in this study with those in other studies 
(e.g., Chand et ah, 2009; de Graaf et ah, 2012; Wilcox, 2012). For example, (Chand et ah, 
2009) used CALIOP in combination with Terra -MODIS observation to compute the DRE 
over the South East Atlantic Region. It is known that low-clouds in this region have a 
strong diurnal cycle driven by solar cloud absorption (Bergman and Salby, 1996; 
Rozendaal et ah, 1995; Wood et ah, 2002). As a result the cloud properties observed by 
Terra-MODIS can be significantly different from those observed by Aqua-MODIS in this 
region, which could lead to different DRE even if the same method was used. 
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5. Summary and Discussion 

Recent advances in satellite-based remote sensing, in particular the launch of the 
space-borne lidar CALIOP, have provided an unprecedented opportunity for studying the 
radiative effects of above-cloud aerosol (ACA). However, the methodologies used in 
recent studies for computing the ACA DRE appear to be either oversimplified (e.g., 
Chand et ah, 2009; Oikawa et ah, 2013) or too cumbersome (e.g., Meyer et al., 2013). 
This paper describes a novel method recently developed for computing the shortwave 
DRE of above-cloud aerosols over ocean. Our method has several unique features 
compared to previous methods: 1) It takes sub-grid scale cloud and aerosol variation into 
account in DRE computations, similar to (Meyer et al., 2013); 2) it treats the overlap of 
aerosol and cloud rigorously by utilizing the joint histogram of COD and CTP in the 
MODIS level-3 cloud product; 3) it relies on grid-level cloud statistics (i.e., COD-CTP 
joint histogram), instead of pixel-level products, and utilizes pre-computed look-up tables 
for ACA DRE computations, making it thus much more efficient than pixel-level 
computations. As shown in Figure 3, Figure 4 and Table 2, DRE computed using our 
method agrees well with the pixel-level computations of (Meyer et al., 2013).. 

In addition to the Southeast Atlantic region, we have recently begun investigating 
the DRE of above-cloud light-absorbing aerosols for global ocean. Some preliminary 
results are shown in Figure 5. We first derived the daily grid-level statistics of above- 
cloud AOD and below-cloud COD, as well as the corresponding ACA DRE, using the 
method described above and then aggregated the daily means to annual mean. The 
temporal aggregation is weighted by the number of ACA pixels in each day during 2007- 
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2010. For example, the annual mean ACA DRE in Figure 5c is aggregated from daily 
mean based on the following equation: 


DRE 


1 ^, -{DRE, 


ACA 


ACA 


In 


( 3 ) 


where {DRE.) ^ is the mean instantaneous ACA DRE in each day averaged over ACA 


pixels, N. is the number of ACA pixels in the grid box in each day, and ( DRE ') is the 

annual mean instantaneous ACA DRE shown in Figure 5c. Figure 5a shows a global map 
of the annual mean 550 mn AOD of above cloud smoke and polluted dust derived based 
on 4 years (2007-2010) of CALIOP aerosol and cloud layer products. Similar to 
(Devasthale and Thomas, 2011), we note several “hotspots” of ACA over the Southeast 
Atlantic, the East-Central Atlantic off the western coast of Saharan Africa, the Arabian 
sea, and the North Pacific basin off the coast of eastern Asia. It is interesting to note that 
the ACA AOD over the east-central Atlantic and Arabian Sea is noticeably larger than 
that over the southeast Atlantic and North Pacific basin. Figure 5b shows the annual 
mean below-aerosol COD derived from the MODIS daily level-3 cloud product using the 
method described in section 3. A notable feature in the figure is that the below-aerosol 
COD over the North Pacific basin is significantly larger than that over other ACA regions. 
Figure 5c shows the annual mean shortwave DRE at TOA aggregated from daily values 
due to ACA smoke and polluted dust over the global ocean. It is intriguing to see that the 
DRE of ACA over the North Pacific basin is significantly larger than that over the 
southeast Atlantic, which is in turn larger than the DRE over the east-central Atlantic and 
the Arabian Sea. In fact, some negative DREs are observed in the latter two regions. This 
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probably due to the COD of below-aerosol clouds being too thin (Figure 5b) over these 
regions to have significant radiative effect, so that the radiative effect of ACA is close to 
that of clear skies (i.e., negative). This is interesting because the above-cloud AOD over 
these regions is actually larger, while the below-aerosol COD over these regions is 
smaller, compared to their counterparts over the southeast Atlantic and North Pacific 
basin. Therefore, the preliminary results seem to suggest that the variability of DRE of 
ACA is modulated by COD, rather than AOD, although it should be noted that we have 
focused only on the light-absorbing aerosols, i.e., smoke and polluted dust, and assumed 
the same narrowband scattering properties for them as in (Meyer et ah, 2013). Further 
research is needed to study the impact of aerosol type and scattering properties on the 
temporal-spatial variation of DRE on global scale. Nevertheless, the preliminary results 
shown Figure 5 clearly demonstrate the usefulness of our new method for global studies. 

It should be noted that this study, and previous ones using CALIOP observations 
(e.g., Chand et ah, 2008; Meyer et ah, 2013; Oikawa et ah, 2013), are limited by the 
capabilities of CALIOP. Arguably, some aerosols exist above every cloud. However, not 
all ACA can be detected by CALIOP due to its inherent limitations. Some ACAs are 
simply too optically thin to be detected, though their radiative effects are also expected to 
be small. Other situations may also be possible. For example, a confined aerosol layer has 
larger volume backscatter than a vertically stretched layer, even if the total aerosol 
amounts are the same, and therefore is more easily detected by CALIOP. Passive sensors, 
on the other hand, are less affected by the vertical distribution of ACA because they 
observe column-integrated scattering by aerosols. Recently, several novel techniques 
have been developed to detect and retrieve ACA properties using passive sensors. 
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(Waquet et al., 2009) developed a method based on multi-angular polarization 
measurements from POLDER (Polarization and Directionality of the Earth Reflectances) 
to retrieve the AOD of above-cloud smoke. This method has recently been extended to 
include both smoke and dust aerosols (Waquet et al., 2013a). Most recently, (Jethva et al., 
2013) demonstrated the ability of a color ratio method to retrieve the above-cloud AOD 
based on MODIS multiple spectral cloud reflectance measurements. A review of the 
emerging satellite -based observations of above-cloud aerosols can be found in (Yu and 
Zhang, 2013). The capabilities and limitations of the passive techniques need to be 
systematically studied through inter-comparisons and comparison with CALIOP 
observations, but they may provide a complementary perspective on ACA. Recall that 
passive imagers have much larger spatial coverage than CALIOP, which makes brute 
force calculations of the DRE at the pixel level computationally expensive. In this regard, 
our method satisfies the need for efficiency of ACA DRE computations based on passive 
imager retrievals. 

As a final remark, we would like to point out that the ACA DRE discussed in this 
study is still a few steps away from the all-sky aerosol radiative effect ( (DRE} all _ sky ). For 

a given grid box, the (DRE} all _ sk>} can be decomposed into the sum of clear-sky and 
cloudy-sky DRE: 

(DRE) = (1 ■ - f c ) (DRE) + f c ■ f ACA ■ (DRE) ACA , (4) 

where f c is the cloud fraction, l^DRE^j dear is the DRE averaged over the clear-sky 
portion of the grid box, f ACA is the fraction of cloudy pixels with ACA detected by 
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CALIOP or other sensors, and (DRE^J aca is the DRE averaged over all ACA containing 

pixels. It is important to note an implicit assumption made in Eq. (4), that is, when a 
distinct ACA layer is not detected, the DRE of ACA is zero. Different sensors (or 
different retrieval algorithms for the same sensor) may have different sensitivities to 
ACA and therefore provide different estimates of f ACA and DRE ‘"'j IC4 . For example, one 
sensor may only be able to retrieve ACA for optically thick clouds. This sensor would 
retrieve a larger (DRE} ACA , but a smaller f ACA , in comparison with another sensor 

capable of retrieving ACA for all clouds. Therefore, when comparing the ACA or all-sky 
DRE estimated based on different instruments or algorithms, it is important to compare 
both the f ACA and [DRE } ACA terms inEq. (4). 
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2 Table 1 Quality control metrics used for screening the CALIOP aerosol layer product. 


Criterion 

CADscore 

<-30 

Horizontalaveraging 

< 20km 

Extinction_QC_532 

0 or 1 

Feature_Optical_Depth_Uncertainty_532 

<-99.5 


3 

4 


1 



1 
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3 Table 2 Regional and seasonal mean values of instantaneous DRE and RFE based on the 

4 pixel-level computation and the new method. 



DRE [W m 2 ] 

RFE [W m 2 AOD' 1 ] 


Bias adjusted 

Bias adjusted 


(unadjusted) 

(unadjusted) 

Pixel computation 

6.63 (5.92) 

55.97 (50.34) 

New Method 

6.39 (5.77) 

53.77 (50.22) 
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CALIOP profile for a grid box 
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MODIS 

COD-CTP joint histogram 



3 Figure 1 A schematic example to illustrate how CALIOP aerosol layer height information 

4 is used in our method to determine the population of liquid-phase clouds below the 

5 aerosol layer in the MODIS COD-CTP joint histogram. 
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3 



2 Figure 2 A schematic illustration of our fast scheme to correct the COD retrieval bias in 

3 the MODIS cloud product due to overlying aerosol contamination. 
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2 Figure 3 a) Seasonal mean (August/September 2007-2011) instantaneous TOA DRE of 

3 above-cloud smoke and polluted dust based on the pixel-level computations from (Meyer 

4 et al., 2013); b) seasonal mean instantaneous TOA aerosol RFE (i.e., DRE per AOD) 

5 from(Meyer et al., 2013); c) same as a), but based on the new method; d) same as b), but 

6 based on the new method. 
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4 Figure 4 Meridional mean DRE and RFE for the region based on the results in Figure 3. 

5 Lines with cross symbol correspond to pixel computations from(Meyer et ah, 2013). 

6 Lines with square symbol correspond to results based on the new method. 
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Figure 5 a) Annual mean AOD (at 550 nm) of above-cloud light-absorbing aerosols (i.e., 
smoke and polluted dust) derived from 4 years (2007-2010) of the CALIOP 5km aerosol 
and cloud layer products, b) Annual mean below-aerosol COD derived from the MODIS 
daily level-3 COD-CTP joint histogram, c) Annual mean instantaneous TOA DRE of 
above-cloud light-absorbing aerosols derived using the new method. 
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